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et 


e experiments at charmonium production threshold are reviewed, it is found that the contri- 


bution of the continuum process via virtual photon has been neglected in almost all the experiments 
and most channels analyzed. It is pointed out that the contribution of the continuum part may 
affect the final results significantly in Y(29) and w(3770) decays, while the interference between the 
continuum amplitude and the resonance amplitude may even affect the J/~ decays as well as the 
(2S) and 7(3770). This leads to the revise of the analysis of strong and electromagnetic amplitude 
in ~(2S) decays, including (25) + VP which is the long lasting puzzle between J/% and (2S) 
decays. For w(3770) physics, a large constructive interference for light hadron modes and destruc- 
tive interference for DD could be responsible for the discrepancy between the larger cross section of 
inclusive hadrons by direct measurement of efe~ — (3770) > hadrons than the DD cross section 
measured using D single-tag and double-tag method. 


PACS numbers: 12.38.Aw, 13.25.Gv, 13.40.Gp, 14.40.Gx 


I. INTRODUCTION 


It is well known that the ete~ experiments have lots 
of advantages in particle physics study: large cross sec- 
tion, small background, and well-determined initial state 
(both four-momentum and quantum numbers). There 
were lots of ete~ experiments, there are still many e+e 
experiments working, and there will be further experi- 
ments build to continue the experimental study. With 
the energy ranges from m*+a~ threshold up to TeV s- 
cale, these experiments contribute a lot to the knowl- 
edge of the world around us. Among them, there are 
some working at the 7-charm energy region, where the 
JPC = 177 charmonium states are produced and stud- 
ied, including Mark-I, Mark-IT, Mark-III, DM2, Crystal 
Ball, BES and so on. Recently, CLEO working at CESR 
decided lower its energy from the BB threshold to the 
charm threshold [1], and BES working at BEPC decided 
upgrade both the accelerator and the detector to make 
a factory-like experiment [2], these two experiments will 
reach an extremely high precision in the study of the 
charmed mesons and the dynamics of the charmonium 
states decay in this energy range. 


J/w, the first vector charmonium state discovered in 
1974 [3], gains lots of attention due to its surprising nar- 
row width and strong coupling to ete” state. Since 
then, it has been used as an ideal laboratory for light 
hadron spectroscopy and charmonium decay dynamics 
study, which are the essential tasks of the low energy 
QCD. The attempt of figuring out the strong decays of 
J/w via three-gluon and electromagnetic decays via one- 
photon annihilation reveals the relative phase between 
these two amplitudes may be large, with the help of some 
pure electromagnetic decay modes of J/w like wx® [4]. 
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This is an important information since it indicates there 
would be no interference between these two amplitudes. 
This situation will be further studied in this paper. 


w(2S), the radially excited spin triplet state of J/y, 
has also the narrow nature and strong coupling to et e7 
state, but most impressing feature found in the w(25S) 
study is the abnormal suppression of some decay modes 
compared with the corresponding J/ decays based on 
perturbative QCD predictions. This suppression was 
first observed by the Mark-II experiment in vector pseu- 
doscalar (VP) decay modes like pr and K*K [5], and 
confirmed by BES [6] (referred as “pz puzzle” in liter- 
atures). Moreover, BES also observed the suppression 
in vector tensor (VT) decays of (29) [7]. This has 
led to substantial theoretical efforts in solving the prob- 
lem [8, 9], unfortunately, most of the models were ruled 
out by the experiments, while some others need further 
experimental data to test. 


The ~(3770), currently regarded as the D wave char- 
monium state, lies above the DD threshold, as a conse- 
quence, the OZI allowed decays of w(3770) + DD would 
dominate its decays. This picture has been considered 
true for pretty long period of time, until recently, a care- 
ful study of the old analyses indicates the DD cross sec- 
tion may be lower than inclusive hadron cross section of 
w(3770), or in other words, there are substantial non-DD 
decays of ~(3770) state [10]. 


The three topics in charmonium studies (relative phase 
between strong and electromagnetic decays, “pr puzzle”, 
and non-DD decays of 7(3770)) play important roles in 
understanding the charmonium decay dynamics. In fol- 
lowing sections, we will examine carefully the experimen- 
tal observables and theoretical expectations in charmoni- 
um study in e*e~ experiments, to provide a possibility of 
investigating these problems in a self-consistent way by 
considering the unavoidable background process in e* e~ 
experiment, namely, the continuum process. We argue 
that, for any exclusive decay final states of these charmo- 


nia decay, in some cases, the contribution of this process 
to the amplitude may be very important, while in some 
other cases, although the direct contribution is relatively 
small, the interference between this term and other dom- 
inant amplitudes may contribute a non-negligible part, 
which maybe provide a guideline to solve the three exist- 
ing problems in the charmonium decays. 

The purpose of our paper focuses on locating the nat- 
ural source of the existing problems rather than to offer 
a detailed solutions of them. So we begin with the inclu- 
sive hadronic process to exhibit the experimental effect 
on theoretical cross section, then simple assumption is 
often adopted to estimate the function on exclusive pro- 
cesses from the continuum contribution. At the last of 
the paper, we also studied the experimental condition de- 
pendence of the results in case the interference was not 
considered, which is true for most of the existing experi- 
mental results on J/w, w(25') and w(3770) decays. 


II. EXPERIMENTALLY OBSERVED CROSS 
SECTIONS 


We know that J/y and (25) decay into light hadron- 
s through two interactions: the three-gluon strong in- 
teraction and the one-photon electromagnetic interac- 
tion. There is in general a relative phase between these 
two amplitudes. This is also true for (3770) in it- 
s OZI suppressed decay into light hadrons. These two 
amplitudes and the phase between them are extracted 
from experimental data by several authors for J/~ and 
(2S) [4, 9, 11]. If we denote the amplitude of three- 
gluon by a3, and one-photon by ay, both of them can be 
complex, the decay rate 


a x |asg +a? . (1) 


In ete” colliding beam experiments, the charmonium 
are produced from ete~ annihilation, there is inevitable 
another amplitude 


ete” > 7* —> hadrons (2) 
accompanied with the production of the resonances. This 
amplitude does not go through the resonance, but in 
general it could produce the same final hadronic states 
as charmonium decays do. So there are three Feynman 
diagrams corresponding to the experimentally measured 
cross sections, i.e. the three-gluon decays, the one-photon 
decays, and the one-photon continuum process, as illus- 
trated in Fig. 1. The former two amplitudes are associ- 
ated with the resonance, while the last one is a slowly 
varying function of C.M. energy (ys). To analyze the 
experimental results, we must take into account three 
amplitudes and two phases. Taking the amplitude of 
one-photon continuum as ac, the experimentally observed 
cross section 


o' x |a3gt+a,ta,|? . (3) 


In this paper, for simplicity, we define 
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e F(s) ig 
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where e is the electromagnetic coupling constant, F(s) is 
the form factor, and 


n Cogg 
39 = M2? +iMT ’ (5) 
Ce? 
"= Wedge ’ (6) 


where Cggg and Cy are taken to be real, and M and I are 
the mass and the width of the resonance. We shall use 
this form for general discussions and numerical calcula- 
tions. 
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FIG. 1: Three Feynman diagrams of hadron process. 
The experimentally observed cross sections in e*e7 
collision are modified by the initial state radiation. For 
the narrow resonances J/w and y(285), the observed cross 
sections are also distorted by the energy spread of the col- 
lider. The radiative corrected cross section is expressed 
as [12] 


Lm 
z o s(l— zx 
anels) = | deF(a,s) eme, (7) 
0 
where Born is the Born order cross section. In the up- 
per limit of the integration £m = 1 — 8/8, y/Sm is the 
experimentally required minimum invariant mass of the 
final state f after losing energy to multi-photon emission; 
F(x, s) has been calculated in several references [12-14] 
and II(s(1 — x)) is the vacuum polarization factor. 
The ete~ colliders have finite energy resolution. The 
energy resolution function G(./s, Vs’) is usually de- 
scribed by a Gaussian distribution : 


G(V/s, V's!) = 
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where A, a function of the beam energy, is the C.M. 
energy spread of the accelerator. So the experimentally 
measured resonance cross section, Jexp, is the radiative 
corrected cross section r.c. folded with the energy reso- 
lution function: 


Texp(V/s) = f ona (VGV, Vs)dv/s! : (9) 
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Through out this paper, all the physical quantities to 
be discussed are experimentally observed ones and the 
radiative correction as well as collider energy spread are 
implicitly taken into account. 

In principle, any experiment measures Terp should sub- 
tract the contribution of the continuum part to get the 
physical quantity Cexp, where Ony and Cezp indicate the 
experimental cross sections calculated from Eq.(9) with 
the substitution of o’ and o from Eq.(3) and Eq.(1) re- 
spectively to og in Eq.(7). Unfortunately, up to now, 
most of the experiments just neglect this contribution 
and Cesp = Oz, is assumed for almost all the channels 
studied, at least at J/w and (25). The noticeable ex- 
ceptions are Pee, Tuu and the total width. These quan- 
tities are measured, together with the resonance mass, 
by scanning around the peak of the resonance and then 
fitting the measured curves with the theoretical cross sec- 
tions. In the fitting, the theoretical cross sections always 
include a continuum term[15]. 

The difference between ez, and Ca implies a plau- 
sible paraphrase in high energy physics literatures. On 
one hand, the theoretical analyses are based on Gerp, on 
the other hand, the experiments actually measure 07,,,. 
However, even in the case that the continuum ampli- 
tude is relatively small, such as in u(2S), certain values 
of the phase possibly lead to non-negligible interference. 
For (3770) scan experiment, the inclusive continuum 
hadron cross section is larger than the resonance peak, 
possible interference may contribute a substantial part of 
the observed cross section. 

We now display the effect from the continuum am- 
plitude and corresponding phase for J/w, w(2S') and 
(3770). To do this, we calculate the ratio 


/ 
Oexrp — Terp 


/ 
exp 


k= (10) 


on 


In order to see the effect of the relative phase, the 
magnitude of a3,, a, and a, are treated as input. In 
principle, a3,, a, and ac are different for different ex- 
clusive modes both in absolute value and in the rela- 
tive strength. For illustrative purpose, following assump- 
tion is used for any of the exclusive mode: the squared 
moduli of a3, and a. are proportional to their branch- 
ing ratios of inclusive hadrons B(R — ggg — hadron) 
and B(R — 7* — hadron) given by PDG [16], and the 
squared modulus of a, is assumed to be proportional to 
the Born order ~*~ cross section multiplied by Rhaa 


TABLE I: Amplitude estimation for three charmonium s- 
tates [16]. o is total cross section for resonance R, (R = 


J/w, W(28), and o(3770)). 


J/ Y(25) (3770) 
lasal” 60% ozy 15% oyes) ~ 1% o(a770) 
Jaz] 17% OJJ 2.9% Ow(2S) 3 x 10 OwW(3770) 
la-|? || 20 nb 15 nb 13 nb 


which indicates the hadronic cross section of the contin- 
uum process, and is estimated by pQCD [16]. Table I 
lists these inputs for J/wW, w(2S) and (3770). 


If we use the BEPC energy spread listed by PDG [16], 
OJy % 3100 nb, cy(25) = 700 nb, and oy3770) ~ 8 nb 
are got. Combining with equation (1), (3) and (10), we 
could obtain k as a function of ¢@ and ¢’, whose variation 
is shown in Fig. 2. It can be seen that for certain values of 
the two phases, k could deviates from 0, or equivalently 
the ratio o2,,/Cexp deviates from 1, which implies that 
the continuum process may produce non-negligible effect 
in experimental measurement. 


(a) k for Tip _ 


(b) k for (2S) 


FIG. 2: k as a function of ¢ and ¢’, (a) for J/y and (b) for 
(28). 


III. a. IN ete” > cc PROCESSES 


A. (3770) 


The ~(3770) is a wide resonance with Fot = (23.6 + 
2.5) MeV [16]. The collider energy spread can be ne- 
glected as long as the standard deviation of the ener- 
gy resolution function A is comparatively small, such 
as A < 2 MeV. Its maximum cross section of inclu- 
sive hadrons in ete~ collision is 8 nb, while the con- 
tinuum cross section is about 13 nb. The resonance pre- 
dominantly decays into DD, while the continuum cross 
section mainly goes to light hadrons. The decay rate 
of w(3770) to light hadrons via three-gluon annihilation 
even though small in the total width, e.g. at the order 
of 1%, or partial width of about 230 keV, which means 
that |a3,| = 0.08|a_| (Y(3770) decays to light hadrons via 
one-photon annihilation is three orders of magnitude low- 
er than the three-gluon annihilation given above [16], and 
could be neglected), if it has interference with continu- 
um amplitude, it could bring an interference of maximum 
1.9 nb in the observed cross section to light hadrons. An- 
other interference between the tail of Y(25) + agg and ac 
at /8 = My 3770) brings an interference of maximum 1 n- 
b, while the interference between w(2S) and (3770) at 
Vs = M3770) is small, with 0.07 nb at maximum. The 
continuum cross section of et'e~ — DD, in the simple 
naive quark model, is estimated by 
4 2P 
aun (11) 


with Pp being the momentum of the D or D. After taken 
into account of radiative correction, it is 0.019nb. But the 
interference between a, and ap is 0.79nb at maximum. 
Here ap denotes the OZI allowed strong interaction am- 
plitude which is responsible for w(3770) + DD decays. 

A possible large constructive interference for light 
hadrons and, at the same time, a large destructive inter- 
ference for DD, could be responsible for the larger cross 
section of inclusive hadrons by direct measurement of 
ete — (3770) —> hadrons [17] than the DD cross sec- 
tion measured by Mark III Collaboration using D single- 
tag and double-tag [18]. 

As to the exclusive decays, it could make some of the 
decay modes with small branching ratios more observ- 
able at the resonance. For example, if the missing decay 
modes of (2S) like pr do appear in w(3770) decays, with 
an enhancement factor [19], their on-resonance cross sec- 
tion could be substantially larger than off-resonance in 
ete” experiment. Quantitatively, if B(~(3770) > pr) ~ 
4 x 1074 (or equivalently, oy(3770)+pr * 0.003 nb) as 
suggested in Ref. [19], and o(ete~ — pr) ~ 0.014 nb at 
Born order by the model of Ref. [20]* , then the max- 
imum interference could be 0.011 nb, much larger than 


cels) = 


[x] The same calculation gives pr cross section at the mass of (25) 


the pure contribution from 7(3770) —> pr. Comparing 
the cross sections on and off w(3770) peak, (3770) > pr 
could be seen through the interference with the continu- 
um amplitude. 


B. (2S) 


As can be seen in Fig. 2 (b), the ratio o7,,,/Cexp could 
deviate from 1 substantially. For each exclusive decay 
channels, k could be different, due to the magnitudes of 
ay and agg, have different coefficients and ac, if estimated 
by form factors, are of different functions of the energy. 
This must be taken into account in the fitting of a,, az, 
and the phase in between. In general, with the interfer- 
ence between a, and the resonance, the maximum height 
of each exclusive channel does not necessarily coincide 
with the maximum height of the inclusive hadrons on 
which data are taken. We shall take the u™ u~ channel 
as an example in the next section. 

In (2S) final state analyses, it is noticeable that the 
observed cross sections of some electromagnetic process- 
es, such as (2S) 3 ntr, W(2S) —> wr, and the fa- 
mous puzzling process ~(2') — pr, are three to four 
orders of magnitude smaller than the inclusive hadron 
cross section of the continuum process, which is about 
15 nb. Form factor estimation [21] gives these cross sec- 
tions comparable to the magnitudes off the resonance[22]. 
It implies that a substantial part of the experimentally 
measured cross section could comes from the continuum 
amplitude a, instead of the 7(25') decays, and interfer- 
ence between these two amplitudes may even affect the 
measured quantities further. Therefore it is essential to 
know the production rate of rr, wr? and pr due to the 
continuum process in order to get their correct branching 
ratios of the (25) decays. 

In order to know whether the observed suppression of 
VP and VT modes in 7(2S) decays are due to the ab- 
sence of strong interaction amplitude, or the destructive 
interference between the electromagnetic and the strong 
amplitudes, or just an incidental destructive interference 
between these two and the continuum process in partic- 
ular experiment, the amplitude a, must be taken into 
account. 


C. J 


From Fig. 2 (a), it is seen that the interference between 
the amplitude a, and the resonance is at the order of a 
few percent which is much smaller than that of (25). 


to be 0.015 nb, which is just below the current upper limit of 
the branching ratio 2.8 x 1075 or upper limit of the cross section 
0.02 nb by BES [6]. 


It is also smaller than the statistical and systematic un- 
certainties of current measurements. Nevertheless, for 
future high precision measurements such as the proposed 
CLEO-c [1] and BES-III [2], when the accuracy goes to 
a few per mille level, it should be taken into account. 


IV. THE DEPENDENCE ON EXPERIMENTAL 
CONDITIONS 


In this section, we discuss the dependence of the ob- 
served cross section in et e~ collision on the experimental 
conditions. The most crucial experiment conditions are 
the accelerator energy spread and the beam energy set- 
ting. The former will smear the intrinsic width of the 
resonance so that change the relative contribution be- 
tween the resonance and the continuum, while the latter 
will affect the relative contribution as well as the abso- 
lute correction to the total rate due to the interference. 
The invariant mass cut or equivalent requirement in da- 
ta analysis will also affect the relative contribution of 
resonance and continuum due to the different energy de- 
pendence of the cross sections in radiative correction. 


A. Dependence on collider energy resolution 


For narrow resonance like J/7 and 7(2S), with intrin- 
sic widths much narrower than the energy resolution of 
the current ete~ colliders, we don’t observe their origi- 
nal resonance curve. Instead, what we actually measure 
is the resonance smeared by the finite energy resolution 
of the collider. In Fig. 3, three cross sections are depict- 
ed: the Breit-Wigner cross section, the cross section af- 
ter radiative correction (Eq.(7)), and the experimentally 
measured cross section (Eq.(9)). 

In actual experiments, data are naturally taken at the 
energy which yields the maximum height of the inclusive 
hadrons. This energy is not the nominal mass of the 
resonance but somewhat higher, neither does it coincide 
with the maximum height of each exclusive channels due 
to the interference effect with ae. For comparison, Fig. 4, 
depicts the observed cross sections of inclusive hadrons 
and yt pairs. Two arrows in the figure denote the 
different positions of the maximum heights of the cross 
sections. It is well know that the radiative correction 
reduces the height of the peak and shift the maximum 
height of the resonance peak upwards, and the energy 
resolution of the collider both reduces the height of the 
peak and shifts it more profoundly. Such shift, depends 
on the energy resolution of the collider, in general could 
be different for inclusive hadrons and for each exclusive 
channels. For example, the peak of u* u” curve is shifted 
more than that of the inclusive hadrons, to 0.81 MeV 
above the w(2.S) nominal mass for BEPC energy spread, 
which is 1.3 MeV. 

It is clear that the shape of the observed cross section is 
much different from that of Breit-Wigner. However, the 
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FIG. 3: Total cross section at ~(25): ogw for Breit-Wigner 
cross section, Or.c. the cross section with radiative correction, 
and Cegzp the measured cross section on a collider with A = 1.3 
MeV. 
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FIG. 4: Cross sections in the vicinity of Y(25) for hadron (a) 
and u™ u` (b) final states. The solid line with arrow indicates 
the peak position and the dashed line with arrow the shift of 
the other peak position. 


energy smear hardly affects the continuum part of the 
cross section. So in the observed cross section, what pro- 
portion comes from the contribution of continuum and 
interference is sensitive to the energy spread. The larger 
the collider energy spread is, the more share the contin- 
uum part contributes in the observed cross section. 


B. Dependence on beam energy 


For demonstration of such dependence, we show the 
curve of ut u` channel, for its dynamics is clear and there 
is no unknown parameter. It is similar to those hadronic 
channels in w(2S') decays which only go through electro- 
magnetic interaction, such as w7? and mtr. Since this 
is an exclusive channel, there is interference between the 
continuum and the w(2S) amplitudes. Such interference 
can be seen clearly from the scan of the 7(2S'), as shown 
in Fig. 5. 
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FIG. 5: Contributions of three parts to the cross section: 
dashed line for QED continuum (øf); dotted line for reso- 
nance (o”); dash dotted line for interference(o’); solid line 
for total cross section(o?). 


The relative contribution of the resonance and the con- 
tinuum changes rapidly as the energy changes, and the in- 
terference term between these two amplitudes also varies 
with energy. The latter could change from negative to 
positive as the energy passes across the nominal mass of 
the resonance. In the actual experimental situation, it is 
important to know the beam energy precisely, which is 
needed in the delicate task to subtract the contribution 
of ac. Possible uncertainty and drift of the beam energy 
need to be taken into account in the determination of the 
systematic errors. 


C. Dependence on invariant mass requirement 


The magnitude of the continuum part of the cross sec- 
tion is sensitive to the upper limit of the integration in 
the calculation of the radiative correction, i.e., it depends 
sensitively on the invariant mass cut in the experimen- 
t (Eq. (7)). This is because that the Born order cross 


section of the continuum term goes up as C.M. ener- 
gy decreases. If the event selection uses a very loose 
cut, then the cross section from continuum part could be 
very large. This is particularly true for exclusive channel- 
s, because QCD predicts the form factors to be powers 
of 1/Q?. While the resonance part is not sensitive to 
the invariant mass cut, since the Breit-Wigner formu- 
la serves as a natural cut by itself. Mathematically, in 
the upper limit of the integration of Eq. (7), as long as 
(1 — sm/s) > T/M, the integration is not sensitive to 
the upper limit, where [ and M are the total width and 
mass of the resonance respectively. 

So in the observed cross section, what proportion 
comes from the contribution of continuum and interfer- 
ence is sensitive to the events selection criteria. Many of 
a time, it is not the invariant mass cut directly applied 
to the data, instead it is affected by many cuts, like the 
momentum cut, kinematic fit, collinearity cut and so on. 
In this case, how much is the contribution of the con- 
tinuum and the interference could only be calculated by 
Monte Carlo simulation. Qualitatively, looser invariant 
mass requirement in event selection would increase the 
share of the continuum part of the contribution. 


It is worth noting here that in principle if a, is not con- 
sidered correctly, different experiments will give different 
results to the same quantity, like the exclusive branching 
ratio of the resonance, due to the dependence on beam 
energy spread, beam energy setting, and invariant mass 
requirement in event selection. This point is especial- 
ly important for the time being, since the beam spreads 
for different accelerators are much different, and events 
selection criteria is very different because of the big back- 
ground in the channels analyzed [23]. 


V. SUMMARY AND PERSPECTIVES 


As we have emphasized in the foregoing discussion, 
the amplitude a,, by itself or through interference with 
the resonance, could contribute significantly to the ob- 
served cross sections in e*e~ experiments on charmo- 
nium physics. Its treatment depends sensitively on the 
experimental details, this has not been fully addressed in 
both e*e~ experiments and theoretical analyses based 
on these results. So far, most of the measurements are 
crude, with large statistical and systematic uncertain- 
ties, so this problem has been outside of the purview of 
concern. Now with large J/w and (2S) samples from 
BES-II [24] and forthcoming high precision experiments 
CLEO-c [1] and BES-III [2], the effect of a. needs to be 
addressed properly. 

To study the continuum contribution, the most promis- 
ing way is to do energy scan for every exclusive mode in 
the vicinity of the resonance, so that both the amplitudes 
and the relative phases could be fit out simultaneous- 
ly. In case this is not practicable, data sample off the 
resonance with comparable integrated luminosity as on 


the resonance should be collected to measure |ae|, which 
could give an estimation of its contribution in the decay 
modes studied. The theoretical analyses based on current 
available ete~ data, particularly on ~(2S) may need to 
be revised correspondingly. 

In fact, another way to free ourselves from the effect of 
the continuum is to analyze the decay product from high- 
er energy experiments. For example, J/y decays could 
be measured from ~(2S) > ntr J/y , (2S) could be 
studied from pp annihilation experiments or from B de- 
cays at the B factories. 
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